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Abstract: 

The quantum system of particles in a double well potential is a widely studied and extremely 

useful example for understanding quantum mechanics. This simple system has recently been 

^ ' used in theoretical proposals and related experiments as a way to make quantum logic gates for 
O ■ 

^ ■ ultra-cold atoms confined in optical lattices (periodic potentials arising from the interaction of 

i atoms with standing waves of laser light). Such quantum gates are the fundamental building 

^ ■ blocks for quantum information processing; in this context the regular array of cold atoms 

CN . in the optical lattice serves as the quantum register. The purpose of this paper is to explain 



X 



how this current research can be understood in terms of well-known principles for systems 
of identical particles, and hence furnish an up-to-date example for the teaching of quantum 
mechanics. 



I. INTRODUCTION 

(N 

Tij- ■ A system with two particles in a double well potential has been studied in recent 

,1,2 



experiments • with ultra-cold atoms. The atoms were confined in optical lattice potentials; 



(N 
O 

^ I these are periodic potentials created by standing waves of laser light. These experiments were 
remarkable as they were able to simultaneously prepare a large number of copies of the dou- 
ble well system, each with a well defined and almost identical quantum state; this quantum 
state could then be manipulated. Furthermore, there has recently been a related theoretical pro- 
posal to make a quantum logic gate^ with neutral atoms held in double well potentials, for the 
purposes of quantum information processing. 

In this paper we explain the basic principles underlying this current research using simple 
quantum mechanics, in order to furnish an up-to-date example for the teaching of quantum 
theory. We shall write out the states of the system explicitly in terms of multi-particle wave 
functions, rather than using the more compact formulation of creation and annihilation opera- 
tions (second quantized notation) as in most of the research papers. In order to establish the 
notation and background, we first consider the system of two particles in a single well — the 
reader who is familiar with this background material could skip to section HVl 



II. TWO PARTICLES IN A POTENTIAL WELL 



The well-known example of two particles in the same potential well illustrates the important 
pedagogical point that the energy levels of quantum systems do not depend on whether or not 
the particles are identical. The occupation of those energy levels, however, is influenced by 
quantum statistics if, and only if, the particles are identical. 

The Schrodinger equation for a system of two particles of mass m moving in one dimension 

is 

( 2mdx\ 2mdx2 ) 

where we take the potential V(x) to be the same for both particles. For this case where the 

particles do not interact, ^ can be written as the product of the wave functions for the individual 

particles 

W(XuX2) = HXl)HX2). (2) 

To provide a simple example we shall consider an infinite square well, i.e. V(x) - for < x < 
L and V(x) = oo elsewhere. Note, however, that the form of the potential makes no diff'erence 
to our subsequent arguments. The single-particle wave functions for an infinite square well are 

= Qr„sm(— ^), (3) 

for n integer, where Qr„ is a normalisation constant; and similarly for i/f(x2). It is convenient to 
express the eigenenergies as En = n^Eg where the ground state energy is given by 

The lowest configuration of the two-particle system has energy Ea = 2Eg, and wave function 

= l/fg(Xi)l/fg(X2). (5) 

Here if/g(xi) is the wave function the ith particle in the ground state. Particles in the first excited 
state have the wave function i/zeixi) with energy eigenvalue E^ (where Eg = AEg for a square 
well). The first excited level of the two-particle system has energy Eg^e = Eg + Eg, and the two 
possible wave functions are 

llfg(Xi)llfe(X2) or l^,(Xi)llfgiX2). (6) 

The degeneracy of the first excited state (and higher lying ones) arises from the symmetry of the 
system under the exchange of particle labels, i.e. the energy of two particles of the same mass 
in the same potential does not depend on which particle is in the excited state. To determine the 
efi'ect of an interaction we therefore need to use degenerate perturbation theory. 



A. The effect of an interaction between the particles 



We now consider the effect of an interaction between the particles so that the Hamiltonian 
becomes 



H = - 



2m dx\ 



2m dx\ 



(7) 



The potential V '{x\ - X2) represents the interaction between the two particles, which we shall 
consider as a perturbation. We only consider potentials that are a function of the distance 
between the two particles \xi - X2I, e.g. the Coulomb potential between charged particles. 

At this stage those who have encountered similar situations before, e.g. in solving the he- 
lium atom, are inclined to anticipate the result and exploit the freedom arising from the de- 
generacy to rewrite the two-particle wave functions as linear combinations the original wave 
functions which are symmetric or antisymmetric under exchange of particle labels. This may 
be prompted by the knowledge that both spatial and spin wavefunctions will be either symmetric 
or antisymmetric for systems of identical particles. However, for now we keep the discussion to 
distinguishable particles, to see the natural emergence of symmetric and antisymmetric states, 
without explicitly invoking the exchange symmetry operator. 

The effect of V'{xi - Xj) on the non-degenerate ground state is given straightforwardly by 
the first-order perturbation: 



(8) 



To evaluate the effect of the perturbation on the degenerate states (of eqn[6l) we write a general 
wave function in the basis of the two states as 



¥(JC1,X2) - ail/g{Xi)lf/e{X2) + blf/e(Xi)ll/g(X2). 



(9) 



This projection of the wave function onto two states is a good approximation when the pertur- 
bation is sufficiently small that is does not mix in significant amplitudes of other states. It is 
convenient to represent these two degenerate basis states as vectors: 



l//g{Xi)lf/e{X2) ^ 



and i//e{xi)i//g(x2) ^ 



^0' 



(10) 



In this truncated basis the Schrodinger equation H^" = reduces to a matrix equation (with 
eigenvalue E): 



Egg + J K 

K Eg f + J 
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(11) 



The eigenvalues are determined by the roots of the equation 

{Eg^e + J-Ef-K^ = 0. 



as 



Here J is the direct integral given by 



(12) 



(13) 



(14) 



and K is the exchange integral (it may appear strange to call this the exchange integral when 
the particles are distinguishable, but it has exactly the same form for identical particles). This 
is given by 

K = J iXiJ iX2 ll/giXi)l//l(Xi)tl/*g(X2)ll/e(X2)V '(xi - X2). (15) 

This is real (the symmetry of the integral for interchange of variables xi <-> X2 implies that 
K* = K). Unlike the direct integral, the exchange integral K cannot be interpreted in terms 
of a simple classical picture; attempts to explain it classically lead to nonsensical concepts, 
e.g. 'exchange forces', which are contrived and completely fictional. (For example, the gross 
structure of the helium atom is determined by electrostatic interactions and nothing else.) 
The two eigenvectors and corresponding spatial wave functions are: 



[l/fg(Xl)ll/e(X2) - lffe(Xl)ll^g(X2)} for Eg^^ + J - K, (16) 



a 













-^[llfg(Xi)llfe(X2) + llfe(Xi)lJfgiX2)} fOT Eg^, + J + K. (17) 



l/^^2 

The eigenstates of the Hamiltonian with an interaction between the distinguishable particles 
are antisymmetric and symmetric combinations of the original wave functions, as in the case 
of indistinguishable particles. To reinforce this point one can consider the case where the two 
particles differ slightly in their mass (or other property aff'ecting their energies), and show that 
this tends to the above result as the masses become equal. 



B. A delta-function interaction 



The interaction between two ultracold neutral atoms can be represented, a good approxima- 
tion, by a delta-function potential 

V\xi- X2) = ad{x2- xi) (18) 



where a is a constant, and a > for a repulsive interaction. For this so-called contact interaction, 
the direct and exchange integrals are the same: 

J = K = ^dx\ipg{xf\ipe{x)\^a (19) 

Since J - K = the interaction does not change the energy of the antisymmetric wave function 
in eqn[T6l which is an obvious consequence of the function being equal to zero when Xi — X2- 
By a similar argument the triplet terms of the helium atom have lower energy than the singlets, 
because the expectation value of the electrostatic repulsion between the two electrons is lower 
for the spatial wave function in which the two parts have opposite sign. 



III. IDENTICAL PARTICLES AND SPIN STATISTICS 



When the two particles are identical their spin has an influence on the occupancy of the 
energy levels. In particular the requirements of overall exchange antisymmetry for fermions, 
and symmetry for bosons, introduces a connection between the spin and spatial wave functions. 
Consequently states with different spin have diff"erent energies, e.g. the singlet and triplet terms 
in helium; it is tempting to think that the energies themselves depend on spin, although they do 
not. The wave functions for two fermions have the form 

— V(^) m(A) \T/(A) m(S) (00^ 

I fermions — i space ^ spin space ^ spin K'^^J 

where (A) denotes antisymmetric with respect to exchange of the particle labels 1 2, and 
(S) denotes symmetric. This subset of all the possible two-particle wave functions has overall 
antisymmetry with respect to exchange of particle labels as required for identical fermions. The 
spin wave functions are written out in Table U (right-hand column); there is a singlet for 
and a triplet for ^fj^^ making a total of 4 functions. A well-known example of this physics is the 
helium atom: for the ground configuration (Is^) there is only a symmetric spatial wave function 
^Pspace = lAis associated with a spin singlet; the first excited configuration (ls2s) is split into two 
terms: a singlet (^S) and a triplet (^S) whose energy separation is twice the exchange integral 
given by 

K = -^ fe'AL('-i)'A;s('-2)^^'A2s(ri)«Ais(r2) (21) 
47r6o J J |r2-ril 

where e^/4neo\r2 - ri| is the Coulomb repulsion between the electron at ri and the electron at 
ra; see Foot^ for further details. For two bosons the wave function is symmetric with respect to 
exchange of the particle labels: 

VU _ vi/(S) \i((S) m(A) \T/(A) 

I bosons — I space ^ spin " space ^ spin" 



a) 



b) 



(/)(x) 



C) 




d) 




<x|L>=^{(|)(x) + ^(x)} 
V2 



<x|R>=^{(/)(x)-^(x)} 
V2 



FIG. 1 : The wave functions of a single particle in a double well potential are illustrated using a 
square well with a high barrier for simplicity, a) The wave function (^(x) of lowest energy has 
even parity, b) The lowest wave function with odd parity 0(x) lies slightly higher in energy 
than (f>, and the energy difference between them increases as the tunneling rate increases, c) 
and d) show the combinations of (p and given in eqns [23] and [241 which represent the particle 

localized in the left and right well respectively. 



This brief review of the physics of two particles in a single well serves as an introduction for 
the treatment of a double system in the following sections. 



IV. DOUBLE WELL POTENTM^L 
A. A single particle in a double well 

We consider double wells that have symmetry with respect to inversion (x -x), so that the 
eigenfunctions of particles in this potential have definite parity (the Hamiltonian of the system 
commutes with the parity operator P, i.e. [H, P] = 0, so these operators have simultaneous 
eigenfunctions). The even and odd parity wave functions for a single particle in a double square 
well potential with a thin barrier are illustrated in Fig. [TJ 

The state whose wave function (0) has even parity lies lower in energy than the wave function 
with odd parity — the barrier has less effect on the odd wave function since has a node at 



X = no matter what the height of the barrier, whereas the energy of the even parity wave 
function decreases as the barrier height decreases. For the simple case of square well potentials, 
when the barrier is reduced from oo to zero, so that the well becomes twice its original length, 
the ground state energy reduces to 1/4 of its original value (see eqnH]). A qualitatively similar 
argument holds for any double well (with reflection symmetry); the solution of the Schrodinger 
equation with no nodes has the lowest eigenenergy. 

We can form linear superpositions of the energy eigenstates 

(x\L) = -^[(f>(x) + ^(x)} (23) 

{x\R) = l^[(f>(x) - ^(x)} (24) 

which correspond to wave functions in which the particle is located predominantly in the left, 
or right, well respectively, as illustrated in Fig. \T\ These are not eigenstates of the system, so 
the wave function of a particle initially prepared in the left well (at t = 0) evolves according to 

Hx,t) = -^U(x)e-'^'"'+^(x)e-'^"''] 
V2 ^ ' 

= J-e-'^"^ U{x) + ^{x)e~'^'] (25) 
V2 ^ ' 

where Q. = (E - E)/ti is the tunneling rate between the wells. When Q.t = Inn, for n integer, the 
wave function is proportional to {x\L), and when 0.t = (2n-l);r the wave function is proportional 
to {x\R) corresponding to the particle having tunneled across into the other well. For negligible 
tunneling between the wells (very high barrier) it is possible to describe the particle as being 
either in the left or right well — this is obviously true when the wells are very far apart so that 
the particle is localised in one of the wells. This is the situation at the start of the experiments 
we shall describe. 



B. Two interacting particles in a double well 

The previous sections have prepared all the mathematical tools we need to understand the 
behavior of two particles in a double well. We choose a potential with a high barrier so that 
initially the system is close to the limit of two separate wells. There is an antisymmetric spatial 
wave function (*Fspace) which can be expressed in terms of the Left/Right wave functions of the 
individual particles as: 

^'%ljxuX2) = {{x,\L){x2\R) - (xmixilL)} . (26) 



This can be rewritten using the ket notation to represent ^Iplce- 



^Sce = { \LR) - \RL) } . (27) 

Similarly, the three possible symmetric wave functions (^P^pLe) 

ILL), \RR), and vpW ^ = -L { |L7?) + |7?L) } . (28) 

Writing the 4 spatial wave functions in this way highlights their similarity with the form of the 
4 spin wave functions for two spin- 1/2 particles; both have one ^P*^^^ and three ^^^\ Combining 
space and spin functions with the restriction of overall antisymmetry, as in eqn [201 we find by 
inspection of Table IJ that there are 6 two-particle wave functions for identical spin- 1/2 particles. 



Now consider the eff"ect of a repulsive contact interaction between the particles. Clearly this 
has very little effect on the energy when the two particles are localized in different wells (with 
small overlap of their wave functions), i.e. for *Fspace and ^PspLe; however, for the wave functions 
containing iffR(xi)iJ/R(x2) (i.e. \RR)) and iAl(-^i)iAl(-^2) (i-e- 1^^)), the interaction raises the energy 
by a positive quantity which we denote as U. The integral for U resembles eqn[T9i 

U = Jdx IM^t a (29) 

and analogously for iI/l(x). Figure [2] and Table II] shows the eigenenergies and associated 
fermionic eigenstates for a double well with a repulsive interaction. 



In the following we shall mainly consider the sub-set of states with the particles in different 
wells. These are separated by the energy gap U from the higher lying states | LL)x^~^ and 
\RR)x^-\ 



C. Quantum gate for two fermions 

In quantum computing notation the state in which both particles are spin-up | TT) can be 
written as |00) and conversely | H) = |11). The state of the system is encoded in terms of the 
four basis states as: 

Y = a|00) + Z7|01) + c|10) + Jill) 

where a,b,c and d are complex coefficients. The state |00) corresponds to a spin-up particle in 
each well for which the antisymmetrised wave function (unnormalised) is 

100 ) CX: {XM)\ Ti ){X2\R)\ T2 ) - {Xl\R)\ Tl ){X2\L)\ T2 ) 



TABLE I: The spatial and spin wave functions for two particles can be those which are 
antisymmetric (A) with respect to exchange of the particle labels 1 2, and others are 
symmetric (S). For a pair of fermions these can be combined to give 6 two-particle wave 
functions of the form 'i'lpice^I'lpL or *i'lpLe*I'l^!i> i-C each of the 6 symmetric functions in the 
bottom row is associated with one of the functions in the top row. This subset (of the 16 
possible product states) has overall antisymmetry with respect to exchange of particle labels as 

required for identical fermions. 



vp(A) 
^ space 

^{|L/?)-|/?L)}^YlpL 


m(A) 
spin 


vp(S) 
^ space 

ILL) 

^{|L7?) + |7?L)} = vpW^e 
\RR} 


vp(S) 

spin 

ltt> 

lU) 



This can be expressed in terms of the energy eigenstates found above as 

100) oc (\LR)-\RL))\n), 

oc ^Seel TT). (30) 

The state |01) corresponds to spin-up in the left well and spin-down on the right for which the 
antisymmetrised wave function (unnormalised) is 

101 ) oc <xi|L)| Ti ){X2\R)\ ii ) - {xm\ ii ){X2\L)\ h ). (31) 

This can be expressed in terms of the energy eigenstates found above as 

101) cx \LR)\n)-\RL)\iV, 

- ( ^Sce + ^Sce ) I Ti ) + ( ^Sce " ^JL ) I it), 

°c ^IpL ( I Ti ) + UT ) ) + KL ( I Ti > - UT ) ) . (32) 



and similarly for |10) and |1 1). Thus we write the quantum computational basis states in terms 
of the 4 eigenstates of two spin- 1/2 particles (which are degenerate for zero tunnelling) as: 



100) 



vj/(-) I 

space I I I ' 



|01) { *^spac&^ ^^spaceA^ 



/^^2 



110) = I^Sc^^"^ 



1 1 1 ) ^space I XX)' 



When the barrier is lowered slowly, ^space and ^^spLe adiabatically evolve into *Pa and re- 
spectively, which are no longer degenerate. Further details are given in Hayes et afi- (see also 
Table HH), but to understand the operation of the gate all we need to know that a phase difference 
accumulates between these states as they have different eigenenergies during the gate opera- 
tion. Thus if the barrier is lowered for a certain period and then raised again (to switch off the 
tunneling), the wave function that starts as 1 01 ) becomes 



^spaceA' } / 



(+) 



(33) 



Wit) = e~'^ I W^i^a''^ + e^^'^Wiliy-' ] I ^^2. 



(34) 



Controlling the process so that the accrued phase shift A(p = n causes this state to evolve into 



-lip 



-lip 



{ *l^spaceA^ 



,(+) , >vp(+) ^(-) 
^ spac&a 



spaC' 



spac 



/V2 
/^^2 



(35) 



The global phase factor e"'' is not important; all other combinations pick up this phase too, 
therefore it does not affect the relative phase. It can be seen that a system initialised in 1 01 ) 
evolves into |10) and vice versa which implements the SWAP operation |01) |10). When 
A(p = In the system cycles back to its initial state. More generally, the states | 01 ) and |10) 
evolve into a coherent superposition, while the |00) and |11) states remain unaffected. The 
SWAP operation can be written as the 4x4 matrix: 

^ 1 0^ 



SWAP - 



10 
10 
.0001^ 
for the basis states in eqn[33l 

It is the so-called square-root-of-SWAP gate that performs a useful entangling operation in 
quantum computation; this operation is represented by the matrix such that 



U 



VSWAP^VSWAP 



SWAP- 



Hence we find that 



VSWAP 



' 1 

i(l+i) i(l-i) 
i(l-i) i(l+i) 



V 















1 



Clearly to implement VSWAP the system undergoes 1 /2 of the phase evolution for SWAP. The 
entangled state obtained by operating on |01) with VSWAP is 



The square-root-of-SWAP gate cannot be broken down into a combination of single qubit gates, 
i.e. gates operating on one qubit at a time (whereas t/swAP can). Single qubit gates would be 
implemented in this system by the rotation of the atomic spin orientation in of one of the wells; 
this may be accomplished relatively simply by optical (Raman) or microwave pulses. The 
square-root-of-SWAP gate, used alongside single-qubit operations, is sufficient to provide a 
universal set of gates for quantum computation^. The implementation of these powerful ideas 
with fermionic atoms in optical lattices is described by Hayes et aP. Similar principles can also 
be applied to a system of two bosons (selecting just two out of the possible spin states to form 
a suitable basis). More details of a gate for bosons are given in the theoretical paper of Vaucher 
et afi and the experimental papers^^. 

V. CONCLUSION 

This paper shows how recent experimental and theoretical research on ultracold atoms in 
double well potentials can be understood in terms of simple quantum mechanics. The experi- 
ments have been carried out with bosonic atoms (the isotope Rb-87) for technical reasons but for 
simplicity the fermionic case has been used to expound the principles. This system described 
by Hayes et aP is particularly intriguing for the way in which influence of the nuclear spin on 
the spatial wave functions is exploited. 

For a more general solution, second quantized notation tends to be used, as in most of the 
papers we cite; this notation becomes more efficient in keeping track of exchange symmetry as 
the number of particles in the system increases. 

There is great current interest in using cold atoms in optical lattices for quantum information 
processing, including the direct quantum simulation of strongly correlated many particle 
systems (such as those in Condensed Matter Physics). The simple approach outlined in this 



f^VswAplOl) = ^ ( 101) + 110) ) + ^ ( 101) - 110) ) . 



paper gives an intuitive way of understanding aspects of these systems. 
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TABLE II: The four lowest energy eigenstates of the two-atom double well system, as shown 
in Fig. [21 Only states a and b are used in the quantum logic gate; these are degenerate in the 
limit of an infinite barrier. (In this limit states c and d are also degenerate, both having energy 
U greater than a and b.) As the barrier height is decreased a and b evolve into states with an 
energy splitting between them and a phase diff"erence accumulates which is an essential part of 

the operation gate described in the text. The localized eigenstates {x\L) and {x\R) are not 
appropriate for a low barrier (high tunneling rate between wells), and one should use products 
of the even and odd parity wave functions and ^ (see eqns.l23]andl24l) — these states are 
included for completeness in the right column. 
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FIG. 2: A schematic diagram of the operation of a quantum gate for two atoms in a double 
well potential. The energies of the four lowest states (see Table HIl) of the 2-atom double well 
system are shown as a function of time as the height of the barrier is lowered and then raised. 
For the initial and final conditions (marked A and C) tunneling is negligible, whereas during 
operation (e.g. at B) the barrier is much lower — indeed the barrier can be reduced to zero so 
that both atoms are in a single well for a certain time before the barrier is raised. The phase 
of the gate is the integral of the shaded region between paths a and b (divided by fi). The 
barrier is altered slowly and smoothly to ensure the adiabatic evolution of the states along the 
given paths (otherwise states a and c mix). The condition for adiabatic evolution dictates that 
the gate must take place over a time significantly longer than h/U. The states have been 



labeled to be consistent with Hayes et afi. 



